How Do Glassy Domains Grow? 
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We construct the equations for the growth kinetics of a structural glass within mode-coupling 
theory, through a non-stationary variant of the 3-density correlator defined in Phys. Rev. Lett. 97, 
195701 (2006). We solve a schematic form of the resulting equations to obtain the coarsening of the 
3-point correlator Xsifitw) as a function of waiting time t w . For a quench into the glass, we find 
that X3 attains a peak value ~ t w 5 at t — t w ~ t w 8 , providing a theoretical basis for the numerical 
observations of Parisi [J. Phys. Chem. B 103, 4128 (1999)] and Kob and Barrat [Phys. Rev. Lett. 
78, 4581 (1997)]. The aging is not "simple": the t w dependence cannot be attributed to an evolving 
effective temperature. 



When a system is quenched below an ordering tran- 
sition, domains of the ordered phase appear and begin 
to grow with characteristic size given by the decay 
length of equal-time correlations of the order parame- 
ter. The corresponding issue for the glass transition has 
been examined numerically 043; using susceptibilities 
and correlation lengths that capture the onset of amor- 
phous freezing 0, U0| , but a quantitative theory of these 
observations has been lacking [f|. Length-scale informa- 
tion similar to that obtained from the 4-density correlator 
and related overlap functions 0, Q has been shown 
to be contained in a certain three-point correlator X3(t), 
whose peak value, and the time at which the peak is at- 
tained, diverge p upon approaching the mode-coupling 
glass transition [lfj. 

In this work we present a theory of the coarsening 
of glassy order, using a non-stationary generalization 
X3(t)tw), whose peak value £l(t w ) is the correlation vol- 
ume as a function of the waiting time t w since the quench. 
We formulate our calculation in the framework of the 
fluctuating hydrodynamics of a dense liquid, and obtain 
results using mode-coup ling theory (MCT) [lOHlJ, m a 
schematic approach [13[ ■ Figures [1] - [5] summarize our 
results. We find that £l(t w ) grows without bound for a 
quench into the MCT glass (Fig. [2]), as t^ 5 , and the re- 
laxation time as i°' 8 (Fig. [5]), in agreement respectively 
with the computer experiments of Parisi and Kob and 
Barrat 0] • As effects beyond MCT cut off the transition, 
the coarsening in experiments, simulations, or a com- 
plete theory will cease at long enough times, but typi- 
cal simulations do not explore these asymptotically long 
time scales and can therefore be compared usefully to our 
MCT coarsening predictions. The three-point function, 
Fi g. Qp shows features incompatible with "simple aging" 
d,[lj] but qualitatively similar to Q . For a quench to 
a distance e from the threshold value on the liquid side, 
fl grows to saturation (Fig. 13]), reaching an equilibrium 
value ~ e _1 , with a relaxation time ~ e -1 - 8 (Figure|l|). 

MCT is an analytically tractable approximation to 
equilibrium liquid-state dynamics that yields a glass tran- 
sition in a homogeneous system, which is why it is so 
widely used despite its shortcomings 0, fll| . In or- 
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FIG. 1: The aging of the two-point function. The correlation 
function C(t,t w ) as a function of t, for various waiting times 
tw shown in the legend. The decay with t becomes progres- 
sively slower with increasing t w . The final parameter values 
are T = 1.0 and A = 2.0. Inset: Scaling t — t w by t r yields 
a data collapse in the a-relaxation regime. Such "simple ag- 
ing", however, is not seen in the three-point function, Fig. 
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der to extend MCT to describe non-stationary states 
such as coarse ning we work with a general field-theoretic 
approach ll|, [T5[, taking care not to use results like 
the Kubo formulae and fluctuation-dissipation relations 
(FDR) , which are justified only in equilibrium treatments 
[16|,ll2| ■ We start with the equations of fluctuating hydro- 
dynamics for the velocity and density fields for an isother- 
mal compressible fluid, extended to large wavenumbcrs 
so as to take into account the modes around the struc- 
ture factor peak In order to obtain an equation for 
the density field alone, we eliminate the velocity while 
retaining momentum conservation but ignoring inertia. 
This yields the dynamical equation 



dS Pk (t) 
dt 



K^pkit) = ^ / ' Vk, q 6p q (t)8 Pk . 



hit), 
(i) 



for the Fourier-transformed density fluctuation Spk(t) 



2 



at wavcvcctor k, with V k , q = k • [qc q + (k — q)cfe_ g ], 
A'i = k B T/S k D L and A 2 = k B T/D L k 2 . Eq. jB can be 
viewed as the no-inertia limit of Eq. (4.1) of [16J. Here 
Dl = (C + 4?y/3)/po is the longitudinal damping, where 
C and 77 are the bare shear and bulk viscosities, k B T is 
Boltzmann's constant times temperature, S k and Ck are 
the equilibrium static structure factor and direct corre- 
lation function respectively and the noise fk{t) obeys 



(Mt)h>(t')) = 



2k B T 



p k (t)5{k + k')6(t-t'). (2) 



From diagrammatic perturbation theory fill [lif , we con- 
struct the equations of motion for the correlation func- 
tion C k (t,t w ) = (5pk{t)8p-k{t w )) and response function 

R k (t,t w ) = (dSp k (t)/df^k(t w )y. 

dCk ^ t tw) = - KxCH&tv) + dsD k (t,s)R k (t w ,s) 



dR k (t, t v 



dsS fc (t, s)C k (s,t w ), 



--5{t-t w )-KxR k {t,t w ) 



ds^k{t,s)Rk(s,t w ), 



(3a) 



(3b) 



with D k {t,t') = {2k B T/D L )p k {t)5{t - t') + M k (t,t'), 
M k (t,t') = (Kl/2) !j% A C q (t,t)C k -SJ) and 
£ fe (M') = Kl^Vl q R q {t,t')Ck- q {t^')- The contribu- 
tion to (f3"a)l from the first term in Dk(t,t') vanishes due 
to causality. Franz and Hertz [l3j obtained schematic 
equations similar to (f3"a|) and (|3b[) for the Amit-Roginsky 
model HI. 

How are the input quantities K\ and Vk.q in equations 
(|3a| and ([3b]) defined for the case of a quench? A compar- 
ison with the treatment of Zaccarclli et al. [171 ] is useful 
here. The Vk, q term in (TTJ) and ^ involves the "resid- 
ual interactions" in [l7|. We define our quench to be an 
abrupt increase in the interaction strength, implying that 
Vk, q should be evaluated at the final parameter value. 
To determine K\, which must now be a time-dependent 
quantity as we are dealing with a non-stationary state, 
we insist, as in [l9j . that for r = (t — t w ) <C t w Eq. 
([3]) obeys time-translation invariance and the FDR. This 
leads, after some algebra, to 

=TR k (0) + Ki [ f Vl q C k - q {t,s)x 

JO Jq 

[^C q (t, s)R k (t, s) + R q (t, s)C k (t, s)]ds. (4) 

In fl7| the term corresponding to Ki enters through the 
equal-time density correlator. The latter being time- 
dependent in a coarsening situation, it is natural that 
K\ should change in time. 

To obtain the equation of motion for the growth ki- 
netics of glassy correlations, we look at the behaviour of 
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FIG. 2: The growth kinetics of glassy correlations. The three- 
point correlation function in an aging structural glassy sys- 
tem, from schematic mode-coupling theory. The peak value 
of Xc(t,tw) grows and shifts to higher t pea k with increasing 
waiting time t w . 



our nonstationary generalization of the three-density cor- 
relation function |9| mentioned above. We introduce in 
the free-energy functional a one-body term e(r)p(r), cou- 
pling the density to an external potential e(r) and lead- 
ing on average to an inhomogeneous shift Sm(r) in the 
mean density field (2fjj . We work in the limit where e(r) 
and hence Sm(r) are uniform, so that the Fourier trans- 
form <5mk has non-zero weight Smo only for wavevector 
k = 0. We will see that this suffices for the purpose 
of extracting the correlation volume. The resulting gen- 
eralized Langevin equation for p, to first order in the 
background density Smo which encodes the effects of the 
field, is [H [H 



96pk{t) ■ KmSp k (t)- k -^^6p k (t) 



dt 



D L 



A 



'-f I V k , g 8p g (t)5p k . q (t) + f k (t). (5) 



Let Ck(t,t w ) and R k (t,t w ) denote the Jmo-dependent 
two-point correlation and response functions implied by 
(|5|). As we are working in a non-stationary state, we must 
define separate 3-point quantities analogous to \3 in 
for C k and R k : Xfc(Mw) = dC k (t,t w ) /dSm \ Smo ^ and 
Xfc(Mw) = dRk(t,t w )/dSm \s mo ->o, with equations of 
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FIG. 3: Evolution of three-point function following a quench 
to a point close to but on the liquid side of the transition. 
Growth to saturation of xo(t,t w ) as a function of t for vari- 
ous t w for A = 0.75, corresponding to a quench into the liquid. 
Inset: The two-point correlator, for the same parameter val- 
ues, shows a relaxation time that grows progressively with 
increasing t w but approaches a finite value. 
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FIG. 4: Three-point function and correlation volume in equi- 
librium. Three-point function and correlation volume in equi- 
librium, (a) xc(t) as a function of t for various e = |A — A c |. 
(b) When Xoif) is scaled with Xc an d time with t peak , we ob- 
tain data collapse for large t. (c) Xc ~ e -10 - (d) Xc ~ tpfak 



motion 
R 



+ K l (t) X %{t,t w ) = jf ds-£ k (t,s) X %(s,t w ) 
+ f dst' k (t,s)R k (s,t w )+S k H (t,t w ), (6a) 



dx£{t,t w ) 

at 



9x *£ tw) +K 1 (t) x ?(t,t w ) = fj dsM k (t,s) X %(t w ,s) 
+ [ dsM k (t,s)R k (t w ,s)+ f dsi: k {t,s)x k {s,t w ) 



+ ds£ k (t,s)C k (s,t w ) + S%(t,t w ), (6b) 
Jo 

where t.' k (t,s) = dE k (t, s)/d6m a \ Smo ^o, and M k (t,s) = 
dM k (t,s)/d5m \s mo ^o, M k (t,s) and 1l k (t,s) are quan- 
tities corresponding to M and £ but evaluated in the 
presence of e(r). The expressions for the source terms 
S k (t,t w ) and S k (t,t w ) are given in the supplementary 
information (SI) and the straightforward but tedious 
derivation of (|6al) and (|6b[) will be presented in a sub- 
sequent paper [201 ]. 

Simplified integral equations keeping track of time- 
dependence alone [23l425| have proved invaluable in ex- 
tracting meaningful results from MCT within a manage- 
able calculation. In this spirit, we suppress dependence 
on wavevector k and write the self-energies in §5§ and © 
as M(t,s) = 2AC 2 (i,s) and T,(t,s) = i\R(t,s)C(t,s), 
yielding equations for the two- and three-point cor- 
relation and response functions which we now denote 
C(t,t w ), R(t,t w ), xc(t,t w ) and xn(t,t w )- We solve the 
resulting schematic versions of (J3J) , (J3J) and ©, whose 
detailed forms are given in the SI section, using the al- 
gorithm developed by Kim and Latz [26, 27 1, to give the 
results quoted at the start of this paper. First, aging 



can clearly be seen in the behaviour of C(t,t w ) in Fig. 
[T] Second and more important is the characteristic non- 
monotone behavior of Xc{t,t w ), and its dependence on 
t w and interaction strength A (Fig. [2]). For a fixed initial 
condition corresponding in our schematic approach to a 
liquid with negligible correlations, we examine in partic- 
ular how Xc(t,t w ) as a function of t changes with t w , 
for values of A corresponding to the liquid and the glass 
phase. Recall that A defines the point to which the sys- 
tem is quenched. For A in the liquid phase but close to 
the transition we find, as expected, that xc(t, t w ) attains 
a peak value Xc a ^ a time t peak , with both Xc an< ^ tpeak 
growing with t w but saturating to finite values as shown 
in Fig. [3J The final peak value of xc{t,t w ) grows as 

(A - Xc)- 1 and xc ~ tpeak ( Fi §- H>- Thcsc final values, 
obtained at t w — > oo are the equilibrium values of the 
corresponding quantities. In the notation of [§| we are 
working at qo — > and our results are consistent with 
theirs in that limit. A more detailed comparison with 
Q or pjj , including an estimate of the correlation length 
requires a calculation of the sensitivity of two-point func- 
tions to a spatially varying potential. 

For a quench into the glassy region, A = 2.0, as shown 
in Fig. [5J Xc g rows without bound. In more detail (Fig. 
[5]), Xc ~ C at a tmie tpeak ~ t h w , with increasing t w . 
The exponents a ~ 0.5 in agreement with simulations 
Q and b ~ 0.8. From the two-point correlator C(t,t w ) 
we find a relaxation time t r ~ i pea fc/4 close to the tran- 
sition. Thus our result implies t r ~ t^j , in agreement 
with the numerical experiment of Ref. [4J ; the relation 
between t r and t pea k remains to be tested. Xc measures 
an effective correlation volume, so that its growth is the 
claimed coarsening of glassy structure, and is consistent 
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FIG. 5: Scaling of order-parameter correlations and relax- 
ation time following a quench into the glass: no simple aging, 
(a) If we scale Xc{t,t w ) by Xa{t w ) and time by t peak , the 
three-point function does not follow a master curve: the be- 
haviour at various t w is not the equilibrium dynamics at an 
evolving effective temperature. This behaviour of \c belies 
the expectation of "simple aging" suggested by the two-point 
function in Fig. [1] (b) The peak height Xc{tw) K t\i 2 . ( c ) 
ipcak and t r both grow as iJJ; 8 . (d) t pca k is proportional to 
t r . (e) If we scale xc(t, t w ) by xc(C = 1/e) and plot it as a 
function of 1 — C(t, t w ), as in Parsaeian et al. no data col- 
lapse is found, tw values are shown in the legend, (f) Scaling 
Xc(t, t w ) by Xc an d plotting it as a function of 1 — C(t,t w ) 
shows data collapse in the regime of a-relaxation. Waiting 
times as in (e). 



with the idea of a growing "domain size" . Regardless of 
the precise values obtained.it is significant that our the- 
ory and the simulations of [2h4[ all find a total structure 
factor Xc growing very sublinearly in time. Our scaling 
laws differ quantitatively (Fig. [S£) from those of 0] , per- 
haps because we measure different quantities. However, 
if we scale Xc(t, t w ) by Xc an d plot them as a function of 
1 — C(t, t w ), data collapse is obtained in the a-relaxation 
regime (Fig. [Sf). We do not claim to understand the 
origin of this scaling or, for that matter, that of Rcf. 

A similar calculation [2(| for the three-point corre- 
lation function for a p-spin spin-glass model with p = 3 
finds again a growing Xc > DU t slower than for the present 
problem. 

We emphasize that the ^-dependent properties we ex- 
tract do not correspond to those of an equilibrium system 
at an evolving A or temperature. Had it been so, scaling 
Xc(t,t w ) by Xc(tw) and time by t peak would have given 



data collapse for all t w as for the equilibrium case (Fig. 
SI). Fig. [5] shows the absence of such collapse even for 
larger t w . It would appear that the 3-point correlator 
is more sensitive to departures from "simple aging" , and 
an interpretation in terms of an evolving effective tem- 
perature, than the two-time correlation function 0, 14 1 
(see inset of Fig. [1]). Perhaps the monotone decay of the 
latter masks such deviations or, more likely, xc carries 
additional, independent information. 

We close by summarising the achievements of this 
work. We have shown that mode-coupling theory 
adapted to describe non-stationary states captures the 
key features of the emergence and coarsening of glassy 
order from a liquid. Through the evolution of a three- 
point function we have shown that the glassy correlation 
volume grows as t^ 5 with waiting time t w , slower than 
domain volumes in conventional coarsening, and the re- 
laxation time of the glass grows as i 1 ^ 8 . These theoretical 
growth laws are supported by simulation studies 0, 0] , 
and the broad features we observe are similar to those in 
In an experimental realization, if the quench is be- 
low the MCT transition but above a putative ideal glass 
transition at, say, the Kauzmann temperature Tr, acti- 
vated processes [29( outside the scope of MCT should cut 
off the growth. Presumably a quench below Tk will giv e 
indefinite growth of a different glassy length scale [3fJ- 
32| with a form not pred icted by MCT. In results to be 
presented separately [20( we find further that an imposed 
shear-rate 7 cuts off aging and coarsening at t w ~ I/7 
in the glassy region and t w = min(i r ,l/7) in the fluid. 
Since the relaxation time goes as t°' 8 , this should imply 
that t r or tp ea k should vary as 7 -0 - 8 . We look forward to 
experimental tests of our results. 
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Zaccarelli for discussions. SKN was supported in part by 
the University Grants Commission and SR by a J.C. Bose 
Fellowship from the Department of Science and Technol- 
ogy 



Supplementary information (SI) 



The schematic form of the equations and the source 
terms 



The schematic form of the equations (3a) and (3b), is 
obtained by following the outline given in the paper. The 
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final forms will be: 
dC(t,t w ) 



dt 



dR{t,t w ) 
dt 



- fx{t)C(t,t w )+2X / C 2 (t,s)R(t w ,s)ds 
Jo 

+ 4A / C(t,s)R(t,s)C(s,t w )ds, 



=S(t-t w )- fi{t)R(t,t w ) 



+ 4A / R{t,s)C(t,s)R(s,t w )ds 



where /i(t) is the schematic version of Ki(t): 
H(t) = T + 6A / C 2 (t, s)R(t, s)di 



(7) 



(8) 



Then the source terms can be written in the form 



Sj?(t, t v 



kEsTc k 

Dr 



Rk(t, t w ) — LO k (t)R k (t, t w ), 



S%{t,t w ) = -B-S*c k (t,t w ) - u k (t)C k (t,t w ) (10) 



Eqs. (7a) and (7b) contain two "source" terms which 
we present in detail here. To proceed, let us define 



*W-g J'lvl, 



Xk- q (t,s)\ -C q (t,s)R k (t,s) 



R q (t,s)C k (t,s) ids 



+C k . q (t, s)\ ^ X g(t, s)R k (t, s) + ^C q (t, s)xk(t, s) 



+ X f{t,s)C k (t,s) + R q (t,s)x%{t,s) 



ds. (9) 



The equations for the three-point correlators are also 
schematicised in a similar way as stated in the paper. 
The final schematic forms of equations (7a) and (7b) will 
be 



9xR ^ tw) + n(t) XR (t, t w ) = 4A jT R(t, s)C(t, s) X r{s, t w )ds + 4A jf R(t, s) X c(t, s)R(s, t w )ds 

+4A jf XR (t, s)C(t, s)R(s, t w )ds + S R (t, t w ) 
dxc ^ tw) + p(t) X c(t, t w ) = Axj'j C(t, s)xc(t, s)R(t w , s)ds + 2A j'j C 2 (t, s) XR (t w ,s)d S 

+4A / C(t,s)R(t,s) Xc (s,t w )ds + AX f X c(t, s)R(t, s)C(s,t w )ds 



+4A / C(t,s)x R (t,s)C(s,t w )ds + S c (t,t w ) 
Jo 



(11) 



(12) 



with the source terms given as S R (t,t w ) = [1 — 
uj(t)]R(t, t w ) and Sc(t,t w ) = [1 — w(t)]C(t, t w ) where 
uj(t), the schematic form of oj k (t), is given as 



tu(t) = 12A / C(t,s)xc{t,s)R(t,s)ds 
Jo 

+ 6A f C 2 (t,s) XR (t,s)ds. 
Jo 



(13) 
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